34* International Conference on High Energy Physics, Philadelphia, 2008 

Combined Analysis of Electric Dipole Moments and Lepton Flavor 
Violating Rare Decays * 

Y. Ayazi and Y. Farzan 

Institute for research in fundamental sciences (IPM), P.O. Box 19395-5531, Tehran, Iran 

In the context of general Minimal Supersymmetric Standard Model (MSSM), new sources for Lepton Flavor Violation 
(LFV) as well as CP-violation appear. We show that in the presence of LEV sources, the electric dipole moment of 
the electron (de) can receive new contributions. In particular, de can receive a significant contribution at one loop 
level from the phase of the trilinear j4-term of the staus, ipA-r- We discuss how we can derive information on <I>A^ by 
combining the information on de with that on the LFV decay modes of the r lepton. We then discuss if this approach 
can bo considered as an alternative to the direct measurement of (f>j^^ at ILC. 



1. INTRODUCTION 



As is well-known, nonzero electric dipole moment of elementary particles would indicate CP-violation. In the 
context of SM, there is an established source of CP- violation which is the famous phase of the CKM matrix. However, 
the contribution of this phase to is smaller than 10"'^* e cm [l| which is too small to be probed in any foreseeable 
future [2] . The phases in the neutrino mass matrix can also contribute to de but their contribution is suppressed by 
fourth power of neutrino mass and is quite negligible: 0(10^^^) e cm fsl]. Thus, detection of a nonzero dg at future 
experiments Q would open a window on new physics. 

Another class of phenomena that can teach us about new physics are Lepton Flavor Violating (LFV) rare decays of 
charged leptons: i.e., fi — )■ ej, t ej and r iij. It is by now established that the violation of lepton flavor takes 
place in the neutrino oscillation phenomenon; however, if the source of LFV is merely the neutrino mass matrix, the 
rate of LFV will be extremely low Q and below the sensitivity of any search in the foreseeable future. Thus, if the 
future searches record a positive signal, it will be an indication for new physics. 

The scale of the new physics might lie at high energies (100 GeV or higher) but we can learn about the properties 
of the new physics by studying the indirect effects on low energy phenomena such as Electric Dipole Moment (EDM) 
and/or LFV rare decay of charged leptons. If there is a way to check what we have learned from the low energy 
phenomena by direct measurements at high energy labs, the results will be more exciting. The former can be 
considered as a guideline for the latter. 

Minimal Supersymmetric Standard Model (MSSM), which is arguably the most popular extension of the SM, 
introduces several sources for CP-violation as well as sources for LFV which can lead to effects exceeding the present 
experimental bounds. The experimental bounds on Br(.^j — > £ij) and the EDM of the elementary particles constrain 
the sources of LFV and CP-violation, respectively. In the context of MSSM with vanishing LFV sources, the bounds 
from the EDMs on the CP-violating phases have been extensively studied in the literature (for an incomplete list 
see Although Ar (trilinear coupling of the staus in the soft potential) is a LF conserving coupling, in the 

presence of LFV, it can affect the properties of leptons of other generations. In particular, in the presence of LFV, 
the phase of Ar can contribute to de at one loop level [8[ . In Q , the bounds on the LFV elements of the trilinear 
A-couplings from the stability of vacuum was overlooked. In this paper, we take into account these bounds and 
demonstrate that at certain parts of the parameter space, these bounds reduce ambiguities in the interpretation of 
results and helps us to derive conclusive bounds. 
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2. SOURCES OF CP-VIOLATION AND LFV IN THE MSSM 

The phenomenology of MSSM is determined by its superpotential and the soft supersymmetry breaking potential. 
The part of the superpotential relevant for this study is 

Wmssm = -YiE^i Li- Hd- Hu ■ Hd (1) 

where Li, and Hd are doublets of chiral superfields respectively associated with doublet {ui en) and the two 
Higgs doublets of the MSSM. In the above formula, e^j is the chiral superfield associated with the right-handed 
charged lepton field em. The index "j" determines the flavor. We have written the superpotential in the mass basis 
of charged leptons [i.e., Yukawa coupling of the charged leptons is taken to be diagonal). At the electroweak scale, 
the part of the soft supersymmetry breaking potential relevant for this study can be written as 

L^ft^^ = - 1/2 (A/iBB + MaW^W^ + H.c.) 

- ({A.YAj + Ay)eX h-Hd + H.c.) - T- {mjJ,,Tj ~ e^/ (mi^)y^ 

- ml^^ Hi Hu - - { Bh Hu ■ Hd + H.c), (2) 

where the "i" and "j" indices determine the flavor and Li consists of {Di e^i). Notice that we have divided the 
trilinear coupling to a flavor diagonal part (AiYiSij) and a LFV part {Aij with An — 0). Terms involving the squarks 
as well as the gluino mass term have to be added to Eq. ^ . The Hcrmiticity of the Lagrangian implies that mf^^^ , 
rnjj^ and the diagonal elements of and are all real. Moreover, without loss of generality we can rephase the 
fields to make M2, Bh and Yi real. In such a basis, the rest of the above parameters can in general be complex and 
can be considered as sources of CP-violation giving contributions to the EDMs. 
After electroweak symmetry breaking, the slepton mass terms can be written as 




where 



{mi)ij = {vnA^ )ij + {ml)i5,j -f cos 2/3(-- -f sin^ 9w)S^J (4) 



{m]f)ij = {rn^Jij + {rnl)^Sij - rriz cos2/3sin^ Ow^ij (5) 

and 

= rriiiAi - tan/?)(5y + Ai,j (Hd) (6) 

in which tan/3 = {Hu) / {Hd) ■ The sources of LFV are {m\)ij, {m'j^)ij and (TOifl)ij — Aij{Hd) with i ^ j. In the 
absence of LFV, at one loop level, each of Aa can contribute to the electric dipole moment of only the corresponding 
charged lepton Eq,. For example, receives a significant contribution from the phase of at one loop level but if 
{m\)eT — (™|f,)eT — Aer — A-^e — 0, the phasc of At cannot induce any contribution to de at one loop level. At the 
two loop level, imaginary Ar can induce a contribution to d^ but the effect is of course loop suppressed [^,0]. When 
we turn on the LFV terms, imaginary A.,- can induce a contribution to d,, at one loop level which can exceed the 
present bound on de by several orders of magnitude. The effect is demonstrated in Fig. [1] For illustrative purposes, 
in this figure, the off-diagonal elements of ra\ and m\ are inserted on the relevant lines as a small perturbation. 
However, to make the analysis, we use the exact formula for de and do not use mass insertion approximation. The 
formulation can be found in the appendix of 

The strong bound on Br(/i 67) [10] implies strong bounds on the e/i elements of rn\, m\ and Aij. Throughout 
this paper we set the e/i elements of these matrices equal to zero. There are also strong bounds on Br(T — > 67) and 
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Figure 1: A neutralino exchange diagram contributing to d^. The photon can attach to any of the Sl, tl, tr or cr propagators. 
The boxes on the left and right sides respectively depict insertion of {m\)eT and (m|j)T-e- The circles indicate insertion of the 
At vertex and the vacuum expectation value of Hd- 



Br(T nj) IjJ but these bounds are about three orders of magnitude less stringent than the bound on Br(/i — > ej). 
Each of the LEV er and /j,t elements can be sizeable (of order of the diagonal elements) without violating the present 
bounds. However if the re and Tfi elements are simultaneously present, both fi and e flavor will be violated and 
Br(/j, ej) can receive a contribution exceeding the present bound on it. To avoid such a situation, we set all the 
fj,T elements equal to zero so the only sources of LEV in the present analysis are the er elements. 



3. NEW CONTRIBUTIONS TO 4 IN THE PRESENCE OF LFV 

In this section, we explore the effects of Ar on de by presenting figures. To draw the figures, the mass spectrum 



corresponding to the a benchmark proposed in 1^ has been chosen. However, the mass spectrum of the staus has 
been allowed to slightly deviate from that at the a benchmark. Notice that at this benchmark, the lightest stau is 
considerably heavier than the lightest neutralino so stau-neutralino coannihilation cannot play any significant role in 
fixing the dark matter relic density. As a result, a slight change of stau parameters will not dramatically affect the 
cosmological predictions. We have checked for robustness of the results and have found that the a benchmark is a 
typical point in the parameter space that demonstrate the overall behavior for most of the parameter space. More 
figures can be found in Q . 

Eig. [2]shows de versus the sine of cjiA^ for = and various values of {m\)eT and {mj^)eT- To draw the dotted 
line marked with (a), {rn\)eT and {mj^)eT are both taken to be large. As seen from the figure, in this case the present 
bound on (depicted by the horizontal line) puts a strong bound on the phase of At-. However, in the case that 
either (m|)eT- or (m|j,)er is very small (as in the case of dashed line (b) and solid line (c)), the bound is considerably 
relaxed. Eigure [T] demonstrates the reason: In order for Im[AT-] to contribute to de, both {m'j^)eT and {m'j^)er have 
to be sizeable. Suppose in the future, rare decay t — > 67 is detected which means "some" of the er elements are 
nonzero. By measuring only Br(r — s- ej), one cannot determine the ratio {rrij^) er / (rnj^) er ■ However, if the number of 

le detected events is statistic 

istribution of the final parti( 
Eollowing 1^, let us define 



the detected events is statistically significant, it will be possible to derive more information by studying the angular 
distribution of the final particles in the r ^ e7 decay [13| . 



, , Jo dcose "^OSC ^^^gg ULUbU 

^P=4x — . (7) 

i (r — > 67) 

where 9 is the angle that the momentum of e makes with the spin of r. In principle, Ap can be measured by studying 



the angular distribution of the final particles at an e e^ collider such as a B- factory 13[. Ap is a, measure of the 
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Figure 2: de versus sinc^A^- The input parameters correspond to the a benchmark proposed in (l^ : = 375 GeV, 
mo = 210 GeV, M1/2 = 285 GeV and tan/3 = 10. We have set |4^|=500 GeV. AU the LFV elements of the slepton mass 
matrix are set to zero except {m'^)eT and (m|j)eT. The dotted (pink) line labeled (a) corresponds to (m|_)eT=3500 GeV^ 



and (m|j)eT=15000 GeV . The dashed (green) line labeled (b) corresponds to (r 



=50 GeV^ and (m|)er=37000 GeV^ 



The solid (red) line labeled (c) corresponds to {mj^)eT=S500 GeV and (mfl.)eT=30 GeV . The horizontal doted line at 



1.4 X 10"^^ e cm 



depicts the present experimental limit [l^ on d. 



hierarchy between left and right LFV elements. That is if {m?j^)^r ^ (m|)er and A^r ^ ^re, converges to —1. 
In the opposite case that {rn^^eT S> {m\)eT and A^t 3> Are, Ap converges to 1. Figs. dSHT]) examine the correlation 
between Ap and d^- To draw these plots we have assigned random values to the er elements of the slepton mass 
matrix in the range satisfying the present bound on Br(r ej) In Figs we have set \At\ = 500 GeV and 

assumed maximal value for the CP-violating phase: (f)A^ = 7r/2. For the LF conserving parameters, we have taken 
the spectrum of the a benchmark 1^. Each pair of the scatter plots shown in Figs. (I3][5|) corresponds to different 
configurations of the er elements. To draw each pair, we have assigned random values (from a logarithmic scale) 
to various er elements. We have then calculated the corresponding values of Br(T 67), \de\ and Ap and have 
depicted the corresponding scatter points with the same color and symbol in fig ures (a) and (b). The horizontal 
lines at 1.4 x 10^-^^ e cm and 10^^^ e cm respectively show the present bound [lO| and the reach of the forthcoming 
experiments y. As seen from the figures, for a given value of Br(T — > 57), de cannot exceed a certain value. 

In the case of Fig. [3l Aer and Are are set equal to zero. As seen from the figure, for a significant portion of 
the parameter space, de lies above the present bound (the points shown with green (light grey) dots). The scatter 
points depicted by pink circle, which appear in Fig. [5]-b as two pink horizontal lines at Ap ~ ±1, correspond to 
de < 10~^^ e cm. From Fig. [3]-b we conclude that for Aer = Are = 0, the bound on de can be satisfied if either 
Br(r 67) is very small (which means that all the LFV masses are very small) or Ap is close to ±1 (meaning 
that there is a hierarchy between the LFV elements). In other words within this scenario, if future searches find 
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Figure 3: a) Scatter plot of de versus Br(r 67). The input parameters correspond to the a benchmark proposed in 'l2{ : 
\^J.\ = 375 GeV, mo = 210 GeV, Mj/a = 285 GeV and tan/3 = 10. We have however set 0a, = 7r/2 and \Ar\ = 500 GeV. All 
the LFV elements of the slepton mass matrix are set to zero except {m\)^r and (mfj)eT which pick up random values at a 
logarithmic scale respectively from (5.9 x 10"" GeV^ 5.9 x 10^ GeV^) and (3.7 x 10"^ GeV^ 3.7 x 10* GeV^). The horizontal 
line at 1.4 x 10~^^ e cm depicts the present experimental limit [l^l and the one at 10~^^ e cm shows the limit that can be 
probed in the near future |2|. b) Scatter plot of Ap versus Br(r — > 67). For each scatter point in Fig. [3]-a there is a counterpart 
in Fig. [3]-b corresponding to the same input values for the er elements which is shown with the same color and symbol. Notice 
that points shown in pink (corresponding to de < 10~^^ e cm) all lie on the horizontal lines at Ap=±l. 



5 X 10 < Br(r e-f) and —0.9 < Ap < 0.9, the bound on de should be interpreted either as a bound on (j)A^ or 
as an indication for a cancelation between different contributions from (jjA^ and other possible CP-violating phases. 

To draw Fig. |4l {m\)(.T and {mj^eT are set equal to zero and instead random values within a range are assigned to 
("^i_R)e-r and {iTij^p)re- The upper limit of the range {i.e., 780 GeV^) saturate the constraints from the Unbounded 



From Below (UFB) consideration 1J|. Notice that these bounds on (™|^^)er and {m^j^re imply a "theoretical" 
bound on Br(r — > 67). The scatter points at the tilted peak with highest de and Br(T 67) correspond to the 
cases that both (rn\fj)eT and {rnj^j^re are close to the upper limit. Notice that a correlation between de, Ap and 
Br(r — » 67) similar to that in the case of Fig. [3] emerges. That is the points marked with green dots (corresponding 
to de >1.4 X 10^^^ e cm), with blue "x" (corresponding to 10^^^ < de < 1.4 x 10^^^ e cm) and with pink circles 
(corresponding to de < 10~^^ e cm) are respectively scattered from right to left. Notice however that in contrast 
to Fig. [3]-b, Fig. |3]-b includes scatter points with —0.9 < Ap < 0.9 and Br(r — > 67) ~ 10~* that satisfy the present 
bound on de (the points marked with "x" in the plot). In Fig. [51 we have repeated the same analysis with the S 



benchmark [l^. In the case of the S benchmark, the constraint from the UFB is so stringent that for all scatter 
points de <2 X 10^^^ e cm and Br(T 67) < 2 x 10^^. 

In Fig. [51 (rn'j^)eT, (™|f)er, {'^'LR)eT and {ni^j^re all take nonzero random values. Fig. [Sl-a contains features of 
both Figs.[3l-a and[4|-a. The significant point is that setting all the er mass elements nonzero, the correlation among 
Ap, de and Br(r — > 67) becomes weaker. That is, unlike Figs. [51-b and [H-b, Fig. [Sib contains points below the 
sensitivity limit of the forthcoming de searches (points depicted with pink circles with de < 10^^^ e cm) for which 
Br(T 67) > 10^^ and —0.9 < Ap < 0.9. This can be explained as follows. At scatter points for which 



[(»Ttifl,)eT, {mDeT < {mlp)re, {mDer] Or [(m|p)e^, (m|)er > {mlp)re, {mp)er], 



(8) 
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Figure 4: Similar to Fig. |3] except that {m\)eT = (n^/fjer ~ and instead {m\ji)eT{= AeriHd)) and (T7^i_R)Te(= Are^Hd)) 
pick up random values at a logarithmic scale from (7.8 x lO"** GeV'^, 7.8 x 10^ GeV^). For each scatter point in Fig. |31-a there 
is a counterpart in Fig. |4]-b corresponding to the same input values for the er elements which is shown with the same color 
and symbol. 




Figure 5: Similar to Fig. |3] except that here in addition to {ra\)sT and {m\)eT, ("^ii^)eT(= Aer^Hd)) and ('m\j{)re(~ Are{Hd)) 
are also allowed to be nonzero. The values of (m|_)eT and (mj^JeT are randomly chosen respectively from (5.9 x 10^'^ GeV"^, 5.9 x 
10'^ GeV^) and (3.7 x 10~^ GeV^,3.7 x 10* GeV^) at a logarithmic scale, (m^^)^^ and [rn^YijTe pick up random values at a 
logarithmic scale from the interval (0.78 GeV^, 780 GeV^). For each scatter point in Fig. [5]-a there is a counterpart in Fig. [5]-b 
corresponding to the same input values for the er elements which is shown with the same color and symbol. 
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Figure 6: a) Scatter plot of de versus Br(r ^ 67). The input parameters correspond to the S benchmark proposed in [12l |: 
\fj,\ = 930 GeV, mo = 500 GeV, M1/2 = 750 GeV and tan/3 = 10. We have however set (j>A^ = n/2 and l^^l = 1800 GeV. 
{m\)eT and (m|j)eT are zero but {m\ji)eT and {m^jijTe pick up random values from (0.18 GeV^, 1.8 x 10^ GeV'^). The 
horizontal line at 1.4 x 10~^^ e cm depicts the present experimental limit lO'] and the one at 10~^® e cm shows the limit that 
can be probed in the near future 0|. b) Scatter plot of versus Br(r — > 67). For each scatter point in Fig. a there is a 
counterpart in Fig. b corresponding to the same input values for the er elements which is shown with the same color and 
symbol. Notice that points shown in pink (corresponding to de < 10~^^ e cm) all lie on the horizontal lines at Ap—±1. 



Al can be of order of Ar which yields —0.9 < Ap < 0.9 but despite sizeable (j>A^, de is still small. Pink circles lying 
in the region Br(r 67) > 10~^ and —0.9 < Ap < 0.9 correspond to such configurations. As a result, without 
independent knowledge of the ratios of LFV elements, we cannot derive any conclusive bound on (j)A^ ■ The fraction 
of the scatter points with de < 10~^® e cm (pink circles) lying in the region with —0.9 < Ap < 0.9 and 10~® < 
Br(r —^ 67) strongly depends on the choice of the range and scale of random pick up of the LFV input. For example, 
had we chosen the lower limit of the range of {m\)eT and {m?p)eT two orders of magnitude higher [i.e., {m\)(.T S 
(0.59 GeV^5900 GeV^) and (m|j)er e (3.7 GeV^ 37000 GeV^) instead of (m|)er € (0.0059 GeV^ 5900 GeV^) 
and (TO^)er G (0.037 GeV^, 37000 GeV^)], no pink circles would have in practice appeared in this region. This is 
understandable because for a constant number of the scatter points, decreasing the lower limit of [m\)f,T and {m'^p)eT 
increases the weight of the region for which condition in Eq. ([5]) is satisfied. 

If by some theoretical consideration we exclude the possibility of conditions (O, the correlation between Ap and 
de is maintained so, for Br(r — > 67) 10~* and —0.9 < Ap < 0.9, the present bound on de can be interpreted as a 
strong bound on (f>A^ For example within the scenario described in 15 1 which relates all the LFV elements to the 
Yukawa couplings, conditions ([8]) cannot be fulfilled. Moreover, in some parts of the parameter space, by combining 
information from different observables with the UFB bounds on the LFV elements of we can exclude the 

possibility of Eq. ([5]). This is demonstrated in Figs. [5] and [T] As seen from Fig. [51 at the 5 benchmark, the bounds 
from UFB exclude the possibility of a contribution from {m\)eT and {nn?p)eT to Br(T 67) at the level of O(10~*). 
Thus, if the 5 benchmark is established at LHC and Br(T 67) turns out to be of order of 10~®, we will conclude 
that the contributions comes from {m\)eT and/or (mj^)er- Moreover if observation shows that —0.9 < Ap < 0.9, 
we will conclude that the contributions of (m|^)e^ and (7^|j)er are comparable so the conditions in Eq. ([8|) cannot 
be fulfilled. In fact. Fig. [7] shows that at the S benchmark, the correlation between Ap and de is maintained even 
when all the er elements pick up nonzero values within the allowed region. Fig. [7] shows that for Br(r — > ej) > 10~^ 
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Figure 7: a) Scatter plot of de versus Br(r ^ 67). The input parameters correspond to the 5 benchmark proposed in 'l2| |: 
l/il = 930 GeV, mo = 500 GeV, M1/2 = 750 GeV and tan/3 = 10. We have however set = 7r/2 and l^^l = 1800 GeV. 
{m\)eT and (m|{)eT respectively pick up random values at a logarithmic scale respectively from (0.23 GeV'^,2.3 x 10^ GeV^) 
and (0.33 GeV^ 3.3 x 10^ GeV^). {mla)er and (mijj)^,; pick up random values from (0.18 GeV^ 1.8 x 10^ GeV^). The 
horizontal line at 1.4 x 10~^^ e cm depicts the present experimental limit [l^ and the one at 10~^^ e cm shows the limit that 
can be probed in the near future b) Scatter plot of versus Br(r — > 67). For each scatter point in Fig. a there is a 
counterpart in Fig. b corresponding to the same input values for the er elements which is shown with the same color and 
symbol. Notice that points shown in pink (corresponding to de < 10~^^ e cm) all lie on the horizontal lines at j4p=±l. 



and —0.5 < Ap < 0.5, (pA^ = 7r/2 yields de higher than the present bound: de > 1.4 x 10~^^ e cm. Moreover 
for Br(r ej) > 10^^ and —0.9 < Ap < 0.9, (j)Ar — yields de detectable in forthcoming experiments: 
de > 10"^^ e cm. 



4. CONCLUSIONS AND OUTLOOK 

We have shown that in the presence of the er LFV elements, the phase of At, (pA^ , can contribute to de at one loop 
level. For values of Br(T — > ej) close to the present experimental bounds, the contribution of (pA^ to de can exceed 
the experimental bound on de by several orders of magnitude. We have found that even if Br(T — > 57) is three orders 
of magnitude below the present bound, the contribution to de can still exceed the present bound on de- The effect 
of (pAr on de strongly depends on the ratios of the LFV slepton masses (jrij^) er / {"iTip) er and {'fnj^p^er / {'fn\p)re- In 
other words, for a given value of Br(T — > 67) and (pA^i Mel can take any value between zero and a maximum which 
depends on the value of Br(T — > 67) and (pA^ [see Figs. (l3]-a)-(l7]-a)]. We have shown that for the specific case that 
(™ifi)eT = ('^i_R)i-e = (see Fig. [3]-b) or (mDer = (™fl)eT = [see Fig. |3]-b], by measuring the asymmetry Ap 
defined in Eq. ([7|) we can solve this ambiguity. However, in the general case that all the er elements are nonzero, as 
shown in Fig.[5l the correlation between Ap and de becomes weaker and to solve the ambiguity, extra information is 
needed. 

Let us suppose that Br(r — > e7) turns out to be close to the present bound {i.e., Br(r e^) > 10^^) and moreover 
let us suppose is measured and found to be —0.9 < Ap < 0.9. Excluding the possibility of a fine tuned cancelation 
between the contributions of different phases, two possibilities emerge: 1) (pA^ is smaller than 0(0.005); 2) (pA^ is 
large but one of the conditions in Eq. ([5]) is fulfilled. To derive a conclusive bound on cpA^. , the second possibility has 
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to be excluded. We show that at some parts of the parameter space such as the S benchmark, the second possibility 
is excluded by bounds on (™|^fl)er and {m\j^T-e from the UFB consideration. 

In summary, combining the information on df. and LFV r decay modes gives invaluable information on . In 
certain parts of the parameter space {e.g., the 5 benchmark), by studying these observables, we can constrain <j)A^ 
however in other parts {e.g., the a benchmark) drawing conclusive bounds on <j)A^ is not possible. In the latter case, 
this method cannot replace the direct measurement of <j)A^ at ILC. On the other hand, direct measurement of (j)A^ 
at ILC can help us to resolve the degeneracies in the pattern of the LFV elements. 
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